Mean field description of the Dicke Model 



Jorge G. Hirsch, Octavio Castafios, Ramon Lopez-Pena, and Eduardo Nahmad-Achar 
Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico 
Apdo. Postal 70-543, Mexico D. F., CP. 04510 

A mean field description of the Dicke model is presented, employing the Holstein-Primakoff re- 
alization of the angular momentum algebra. It is shown that, in the thermodynamic limit, when 
the number of atoms interacting with the photons goes to infinity the energy surface takes a simple 
form, allowing for a direct description of many observables. 
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I. INTRODUCTION 

The Dicke Hamiltonian has been successfully employed to describe the superradiance when two-level atoms interact 
with a one mode electromagnetic field [T] , in superconducting circuit QED systems [3] , and in the open experimental 
setup [2] where the phase transitions in a superfluid gas inside an optical cavity and circuit QED was observed ^ [S] . 
The Holstein-Primakoff realization of the quasi-spin operators representing the collective atomic excitations has been 
widely employed in its description [3HS]- 

The expectation values of the energy per particle, the number of photons per particle and the fraction of excited 
atoms, has been calculated employing a truncated version of the Dicke Hamiltonian in the normal and superradiant 
phases in the thermodynamic limit, and found to be in good agreement with the numerical estimations of these values 
even for a small number of atoms 6J. In this contribution we show that these expressions can be found directly 
from a mean field description of the Dicke model. They also reproduce the expressions found employing SU(2) 
coherent states, but fail in describing the fluctuations in the number of photons and atoms, which requires the use of 
symmetry-adapted coherent states [U [8] . 

The Dicke Hamiltonian describes the collective interaction of N two-level atoms with energy separation equal to 
HQ A with a one mode radiation field of frequency ujp. It is given by 

Hd = a^a + loaJz + ^ (a^ + a) ( J_ + J+) , (1) 
v N 

where uja = loa/C^f > is given in units of the frequency of the field, and 7 = ^ /ojp is the (adimensional) coupling 
parameter. The operators a, a\ denote the one-mode annihilation and creation photon operators, respectively, 
the atomic relative population operator, and J± the atomic transition operators. 

II. MEAN FIELD DESCRIPTION 

In the literature it is customary to employ the Holstein-Primakoff representation of the angular momentum operators 
0, with j = N/2, 

J+ = 6V2j - b'^b, J_ V2j - b^b b, ^ b^b-j, (2) 
where new Bose operators b^,b are introduced, which obey the commutation relation 

[6,6t]^l, 

and the vacuum of the new bosons which satisfies 6|0f,) = is |0;,) — \ j, — j). Making these substitutions into Ho, Eq. 
([1]), the two-mode bosonic Hamiltonian becomes 

Hhp^ aU + UAib^b-j) (3) 

+7(at+a)(6t.yrT|i+yr^fc). 

The mean field description of this Hamiltonian is easily obtained employing as a trial state the direct product of 
two Heinseberg-Weyl coherent states [TU] \a) and which are eigenstates of the two bosonic annihilation operators 



a\a) b\(3) =/3|/3), 



(4) 
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where a and /3 are complex numbers, parameterized a.s a — pa e"^° , P — pb e"^*" , with Pa, Pb > and < (j>a,<t'b < 27r. 
The expUcit form of these coherent states is 



_^^(CaX C 1^^^ with C = «,/3. (5) 



n! 

71=0 



In the above expression \n) is the normahzed Fock state with n bosons. 

To obtain the energy surface £{tL>Aj 7i j] ct, /3) = {a/3\HHp\a/3) some matrix elements are directly obtained from Eq. 
Q. They are 



{a\a^a\a) = \a\^^pl, {/3\b^b\f5) ^ ^ pf, 
{a\a^ + a\a) = a + a* = 2 pa cos(f>a- 



More interesting is the evaluation of the matrix element of the last term in Hhp, which requires the use of the 
explicit form of the coherent state, Eq. ([s]). For the square root to make sense in the last term of Eq. ([3]), the 
maximum number of b bosons must be limited to 2j. 

(/3|6t^l-f + ^1 - f b\(3) = 

2i 



2e-Pt pbcoscj^b E 



V 2i 
71 = V 

Putting together the previous expressions, the final form of the energy surface is 
£{uJA,l:j;a,(3) = pI + (pI - j)ujA+ PaPb cos 0a cos 06 e"''" 22 

71 = » ^ 

III. THE THERMODYNAMIC LIMIT 

From the Holstein-Primakoff representation, Eq. (2]), it follows directly the Heisenberg-Weyl contraction of the 



SU(2) group. This contraction corresponds to the t 
Heisenberg-Weyl coherent states become equivalents [TT] 
The function 



lermodynamic limit j — > oo, and in this case the SU(2) and 



2j 



2n 



has the asymptotic limit 



Lim,^oo^^(p,j) = Jl-^. (10) 



2j 



This expression can be "guessed" doing the substitution 



which is equivalent to replace (/3|(6^6)"|/3) by in the series expansion of the square root. All the terms 

which are not in normal order are neglected in doing so. 

As it can be seen in Fig. jTj where F{p,j) is evaluated as a function of and compared with the simpler limiting 
expression, the convergence is fast, and even for j — 100 it is hard to distinguish between the different calculations. 

In the thermodynamic limit, when the number of atoms goes to infinity, the energy surface £ takes the simple form 



^ ^ pI + {Pb~ j)^A+^lPaPbjl - ^COS0aCOS06 (12) 
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FIG. 1: (Color online) The function F{p,j) evaluated as function of -7^, for j=10 (dashed blue line), 100 (dotted red line) 
and 1000 (green dot-dashed line), and its hrnit for j — )■ cxj (continuous black line) . 



For a given set of Hamiltonian parameters i-oa,J, the values Pai't'acT Pb,4>b which minimize this expression provide 
the mean field wave function, the best approximation to the exact ground state as a product wave function. They are 
obtained by solving the equations for the energy surface critical points 

dE d£ 

7^=0, 7rr=0, t = a,b. (13) 

opi d(t>i 

The explicit form of these equations is 



Po + 27^6^1 - cos 0a cos (/)b = 0, (14) 



7PaP6y 1 - ^sin0aCOS06 = 0, (15) 

I 7 

7PaCOS0QCOS(/)f, - w^PbWl - ^ = 0, (16) 



2j 

And the solutions associated with the minima of the energy surface are 

Pa = 0, (f>a undetermined, 
Ph = 0, (/)(, undetermined, 

0, 



7PaPfc\/l - ^COS (/la sin = 0. (17) 



if 7 < 7c, (18) 
if 7 > 7e, (19) 



Pb 
with 



Substituting these expressions in Eqs. Q and (12), the expectations values of the energy per particle, the number of 
photons per particle and the fraction of excited atoms are obtained, in full agreement with previous works [^{5]. 
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IV. CONCLUSIONS 



We have shown that a simple expression can be found for the energy surface of the Dicke Hamiltonian, employing a 
Holstein-Primakoff realization of the rotor algebra for the atomic sector, and using as trial states the boson coherent 
states associated with the photons and the atoms. A variational treatment of the energy surface allows to obtain 
analytical expressions for some intensive observables, which reproduce those obtained previously with truncated 
versions of the Dicke Hamiltonian [3]. They also reproduce the expressions found employing SU(2) coherent states, 
but fail in describing the fluctuations in the number of photons and atoms, which requires the use of symmetry-adapted 
coherent states [3 [H] . 

This work was partially supported by CONACyT-Mxico, DGAPA- UNAM (projects IN102811 and 102109), and 
FONCICyT (project 94142). 
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